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Reversible compression of an optical piston through Kramers dynamics 
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We study the reversible crossover between stable and bistable phases of an over-damped Brownian 
bead inside an optical piston. The interaction potentials are solved developing a method based on 
Kramers’ theory that exploits the statistical properties of the stochastic motion of the bead. We 
evaluate precisely the energy balance of the crossover. We show that the deformation of the optical 
potentials induced by the compression of the piston is related to a production of heat which measures 
the non-adiabatic character of the crossover. This reveals how specific thermodynamic processes can 
be designed and controlled with a high level of precision by tailoring the optical landscapes of the 
piston. 


Optically trapped Brownian particles constitute ideal 
test systems for non-equilibrium statistical physics with a 
great variety of stochastic protocols under external force 
fields that can be implemented [QQ. A particular atten¬ 
tion has been devoted to measuring thermal fluctuation- 
induced escape over an optical potential barrier and ex¬ 
ploring Kramers rate theory, including the observation of 
stochastic synchronization mg. More recently, quantita¬ 
tive tests of so-called fluctuation theorems have involved 
optically trapped nanoparticles, both in the over- and 
under-damped regimes mni- Bistable optical potentials 
are currently exploited for developing Szilard-types en¬ 
gines and studying the connections between information 
theory and thermodynamics [T014T2] . 

In this Letter, we monitor, at room temperature, the 
crossover between stable and bistable motions of a ther- 
malized over-damped Brownian particle optically trapped 
in front of a mirror. For specific positions of the mir¬ 
ror, the coherent superposition of the incident trapping 
beam and the reflected beam induces dynamical bistabil¬ 
ity where the particle is activated between two distinct 
positions along the optical axis. We demonstrate that the 
whole interaction potential can be solved by interpreting 
the two positions as distinguishable metastable states. 
Diffusion limited escape rates and associated activation 
energies are extracted, together with the actual distance 
separating the metastable states. Remarkably, this is 
performed without any preliminary spatial calibration of 
our optical setup. 

While the instantaneous position of the particle is a 
stochastic process, the position of the mirror is an exter¬ 
nal variable that controls the optical force field applied 
to the particule. We show that the movable mirror acts 
as an optical piston that quasi-statically injects reversible 
work into the system in the form of Helmholtz free energy. 
From this description, the energy cost on the Brownian 
particle associated with the compression of the piston can 
be measured precisely. This reveals that the crossover is 
close to the adiabatic limit. We show that the deforma¬ 
tion of the optical potentials through the displacement 
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of the piston produces some reversible heat that fully 
accounts for this small deviation from perfect adiabaticity. 
Remarkably, this relation between the optical landscapes 
and the production of heat points to an efficient resource 
for designing specific thermodynamic processes. 



FIG. 1. (a) Sketch of the trapping configuration. A micron- 
size bead (radius R — 500 nm) is trapped by a Gaussian 
beam (A = 785 nm) focused through a microscope objective. 
The position of the bead z is defined with respect to the 
beam waist wq. The end-mirror is positioned at a distance £ 
from the waist, (b) Evolution of the optical landscape in the 
vicinity of the waist (z = 0) and as a function of £. Blue colors 
correspond to regions of higher intensity, i.e. deeper potential 
energy and thus to stable positions. The case of a stable 
landscape corresponds to the distance £\ (green line) and a 
bistable landscape is crossed at £2 (blue line). Plots of the force 
diagrams associated with G-stability (c) and G-bistability (d). 


In our experiment, a single polystyrene bead is optically 
trapped by a focused Gaussian beam in a water cell at a 
typical 2 /am distance from a metallic mirror (see details 
in Appendix |A|. The trapping beam, characterized by a 
fixed waist wo located at z = 0, propagates in the fluid 
along the z > 0 optical axis with a wave vector +k -see Fig. 
1 (a). It is Mx magnified through the transparent bead 
acting as a lens and reflected with a reflection amplitude 
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r[ A] by the mirror placed at a distance £ from w^. This 
creates a coherent optical landscape at the positi on o f the 
bead z measured from the waist (see Appendix |B 1 | ) 

ITM = IE +k (z) + r[A] • E ~ M \z - 2£)| 2 , (1) 


displayed in Fig. 1 (b) as a function of 2 and £ for a fixed 
value of M. As expected from its interfering nature, the 
optical landscape profile changes with the waist-mirror 
distance L 

The corresponding evolution has direct consequences 
on the dynamics of the trapped bead in the vicinity of the 
waist. Although the approach we develop below is fully 
general, these consequences are most easily described in a 
dipolar approach. Here, the non-absorbing bead is mod¬ 
eled by a real dipolar polarizability a and the optical inter¬ 
action potential is Um{z, £) = —A • a • 1^ (z, £) / ( 2 £qti 2 c) , 
with 77-2 the refractive index of the fluid (see Appendix 


B 2 ). The coupling constant A allows accounting for bead 


size effects, with i « 1 beyond the dipolar limit B3I- 
Within such an approach, the time-averaged conservative 
optical forces acting on the bead directly derive from the 
potential energies F opt (z,£) = —d z UM{z,i). The dipolar 
approach therefore reveals in a straightforward way the 
crucial property that the optical force field is directly 
determined from the optical landscape, for every choice 
of £. 

Because of the coherent nature of Eq. ( 1 ), one dynami¬ 
cal configuration can be selected from the distribution of 
the successive resonant phase conditions in the (z, £) space 
(corresponding to regions of maximal intensity shown in 
Fig. 1 (a)). This roots the analogy with a piston-like 
action exerted by the mirror on the Brownian bead. For 
instance, picking at t\ a resonant phase condition precisely 
on the waist as shown in Fig. 1 (b) leads to restoring 
forces that will maintain the bead in a stable trapped 
position at z = 0 as displayed in Fig. 1 (c). But a mere 
compression of the piston to £2 brings a bistable configu¬ 
ration where the resonant phase evolution induces regions 
of local stability from both sides of the waist separated by 
an unstable point at z = 0 , as seen in the bistable force 
diagram Fig. 1 (d). 

The stochastic trajectory z(t) of the bead evolves in 
every such ^-configuration, modulating the intensity recol¬ 
lected by the objective (see Appendix [Cj. The time-traces 
of these modulations allow us to retrieve the essential 
features of the potentials explored by the bead. For the 
stable ^-configuration of Fig. 1 (c), the time-trace is 
displayed in Fig. 2 (a) and corresponds to a Brownian 
motion performed in an quasi-harmonic potential probed 
by the bead at the bottom of the whole optical potential 
shown in Fig. 2 (b) -see below. The corresponding trap 
stiffness is measured by a power spectral density (PSD) 
analysis [14|. At low Reynolds numbers, this only relies 
on the determination of the roll-off trapping frequency 
and on the knowledge of the fluid friction coefficient 77 . 
Note that we neglect the ^ 20% systematic error on the 
perpendicular viscosity when working at a 2 /im distance 
from the surface. 
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FIG. 2. Experimental intensity time traces respectively as¬ 
sociated with a Brownian motion in a quasi-harmonic well 
(a) and in a bistable potential between two distant spatial 
positions (c). The plots (in units of k^T) in (b) and (d) are 
the potentials associated with the time traces (a) and (c), 
respectively. 


The time-trace displayed in Fig. 2 (c) corresponds to 
the bistable ^-configuration of Fig. 1 (d). The inter- 
mittency of the intensity signal between the two distin¬ 
guishable mean values is the signature of the activation 
of the bead between two metastable positions along the 
optical axis. There are indeed clearly two different time 
scales: a short one associated with Brownian fluctuations, 
and a much longer one on which take place activating 
events from one to the other of these two positions. In 
each metastable states, similar time-traces as those of 
Fig. 2 (a) reveal a quasi-harmonic motion, expected for 
local equilibrium. Therefore, while the bead performs its 
Brownian motion within a local well, it diffuses across the 
potential barrier through rare events thermally assisted 

m- 

The separation of dynamic time scales and the coherent 
nature of the optical landscape that provides a built-in 
spatial reference are two sufficient criteria for applying 
Kramers’ theory to our problem m- In this framework, 
the interaction potential of the bead is reconstructed for 
any length of the optical piston, without resorting to any 
position density probability of the bead along the optical 
axis that would require an absolute spatial calibration of 
the setup. As soon as the process is stationary with a 
sufficient number of recorded activating events, Kramers’ 
theory connects escape rates evaluated from averaged 
residency time T{ within each {i = 1 , 2 } well [T7] 

_ 1 _ _ f UM(zbJ) ~ Um(zu£) \ 

Ti 2 ttt] exp V k bT ) 

to local trap stiffnesses Ki = d^UM(z, £)\ Zi that fix the cur¬ 
vature at the bottom of each well, and to the actual shape 
of the barrier (position 27 , and height) through the absolute 
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£ (jum) 

M 

A 

A z (nm) 

Ubi ( k b T ) 

ftb (pN//rm) 

fti (pN//xm) 

ft 2 (pN/fim) 

AU (kbT) 

n (s) 

r 2 (s) 

1.811 

1.505 

2.92e-3 

285 

2.98 

2.81 

3.35 (3.12) 

3.51 (3.85) 

-7e-4 (-7e-4) 

0.505 (0.526) 

0.494 (0.474) 


TABLE I. Parameters of the interaction potential extracted from the resolution method applied to the bistable .^-configuration. 
Values directly obtained experimentally are indicated in brackets. 


value of its curvature ^5 = —<9 2 Um(^ ^)\z b - Taking the ra¬ 
tio of both rates therefore leads to measuring the potential 
energy difference between the local equilibrium positions 
A U = U M (z 2 J) - U m (zi,£) = k B T In (n/r 2 \[k\ /k 2 ). 

Measured ftp2 and t\^ provide a non-linear system 
of equations which solution fixes the three (M, £, A) 
parameters needed for a de finition of the interaction po¬ 
tential (see Appendix |B 3| ). Experimental values hav¬ 
ing their own uncertainty, the precision on AU is below 
fepT/2 and below 6 nm for £ (see Appendix HI- We also 
extract from the resolution algorithm the barrier posi¬ 
tion Zb, inverted curvature ft& and height, measured as 
Ubi = — Um{z%, tj. The barrier, 3 times higher 

than /cpT, is still shallow enough to allow the bead map¬ 
ping, through thermal fluctuations, the bistable potential 
around z = 0. The distance A z over which the bead is 
activated is also measured. From the parameter values 
gathered in Table I, the interaction potential profile can 
be plotted as a function of the bead displacement as in 
Fig. 2 (d) in units of k^T. We stress that the phase 
structure of 1°^ forbids a simple 4 th -order potential (i.e. 
Duffing type). 

Further insight comes from looking at the system from 
the point of view of the evolution of configurations con¬ 
trolled by the external variable £. An incremental change 
d£ of the length of the piston pushes the bead out of 
equilibrium and forces it to relax in the new £ + d£ con¬ 
figuration with a stiffness ft and a fluid friction 77 on a 
time tD ^ ft/77 set by diffusion, typically ca. 10 -2 s in our 
conditions. Hydrodynamic effects on the bead due to the 
motion of the mirror can be neglected since the incremen¬ 
tal shift of the mirror by \d£\ =20 nm, performed with a 
speed of 1 mm/s set for the piezo-actuator, is associated 
with a low 10 -8 Reynolds number. Accordingly, the dis¬ 
placement of the fluid remains purely diffusive and the 
moving piston therefore has no direct mechanical action 
on the bead. Under such conditions, the only source of 
mechanical loss in the system is given by the relaxation 
process from one configuration to the other. 

We emphasize that z(t) does not map the entire 
canonical equilibrium distribution associated with an £- 
configuration. It only maps a thermally accessible subset 
of it, that can be resolved for sufficiently long acquisition 
times (30 s in our experiment). The notion of stabil¬ 
ity then corresponds to local stable wells much deeper 
than k^T while bistability corresponds to local barrier 
heights of the order of k^T over which the bead can be 
activated. In this picture, stable phases can be identified 
from bistable phases, as drawn in Fig. 3 (a) in the M — £ 
parameter space. 

Our coherent optical piston configuration gives a unique 


capacity in monitoring the crossover between these phases. 
Indeed, a continuous compression of the piston connecting 
two stable configurations forces the bead to go through a 
whole phase of bistability, starting for a piston length £i 
with the bead in an initial stable position at the incident 
waist z = 0 and ending for £f < £{ with the bead in the 
very same spatial position but within a different stable 
potential. From our resolution method, the bistable dy¬ 
namics of the bead can be solved for each step in £. The 
(M, t, A) parameters fixing each potential throughout 
the bistable phase are extracted and the actual steps of 
the entire path followed by the system between £i and £f 
can be plotted in the M — £ plane. It is worth noting that 
the resolved £ values follow precisely the mirror actuation 
command and that the path shows a small dispersion 
in M values. This important outcome of the analysis 
leads to determine the potential profiles even in the stable 
phase from a mere extrapolation on the variable £. This 
is done for instance for the stable £\-configuration of Fig. 
1 with the measured potential profile plotted in Fig. 2 
(b). 

As shown in Fig. 3 (b), the path can also be represented 
through intensity probability densities, associating to 
each probability density extrema a well of the resolved 
potentials. These plots clearly reveal the progressive 
onset of a bistable dynamics of the bead along the optical 
axis of the setup. As a clear advantage of our statistical 
method, this crossover dynamics can be probed despite 
the unknown exact relation between measured intensities 
and bead positions along the optical axis. 

It is possible to give a thermodynamic description of 
the path from an incremental energy balance. Following 
m, this can be drawn from the Langevin over-damped 
dynamics of the bead accounting for the contribution of 
the piston with 

d U m {z,£) = dQ + dnU M [z,£)d£, (3) 

where d Um{z,£) is the change of the potential energy 
of the bead. This change comes from two sources: (i) 
the heat balance dQ = £dz — rjzdz between fluctuation 
(determined by a thermal stochastic force £) and friction 
(related to 77), and (ii) the external work done on the 
bead by the displacement of the piston where the ex¬ 
ternal variable £ controls the evolution of the potential 
configurations. 

By waiting much longer than tp between each incre¬ 
mental change d £, steady-state of every new configuration 
is reached through mechanical relaxation of the bead. 
This insures that the system evolves through the configu¬ 
rations in an iso-thermal and quasi-static way. Moreover, 
all incremental changes in the optical potential are kept 
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smaller than k^,T. This implies that the bead would go 
back exploring the same configurations if the displacement 
of the piston would be reversed. As a consequence, the 
whole path in Fig. 3 is thermodynamically reversible. 

Under such conditions, the averaged external work 
is directly related to the Helmholtz free energy with 
(deUM(z,£)){ z } d^ = d F(i) [Till [20] . The averaging pro¬ 
cess being performed over the positions occupied by the 
bead in the given configuration, the free energy is only 
function of the external variable i. The total amount 
of work performed by the piston through the isothermal 
reversible process is directly given by the free energy dif¬ 
ference between the initial ii and final if positions of the 
piston which can be calculated from the ( ii,if) canonical 
partition functions as 


W = F(i f ) - F(li) = -k B T\n 




( 4 ) 


For both initial and final stable configurations, the par¬ 
tition functions can be expanded to second order around 
the waist as 


^m( 0 ,q,f) 

Z{iij) ~ e 


/ 2tt/c r T 


^ij 


( 5 ) 


where Kij = S 2 Um(£ 5 ii,f)\z~o are the stiffnesses of 
the (ii,if) stable potentials. The total energy balance 




Magnification M 


Intensity (V) 


A U M = W ± Qrev, with 


A Um = CTm(0,^/) — Um( 0, A) 


Q rev k B T In 



(6) 


connects, along the path, the total amount of energy 
change A Um to the heat Q rev produced by the whole 
reversible process. For the crossover of Fig. 3, A Um = 
—2.34 x 10 -21 J (±3%). As clearly seen, a stiffness 
difference between the initial ii and final if configura¬ 
tions is directly related to the production of heat. We 
unambiguously calculate a quantity of reversible heat 
Qrev = —2.48 x 10 -22 J (±10%) transferred to the fluid 
by the bead along the path (see Appendix [E] for the eval¬ 
uation of the uncertainties). The negative Qrev value 
means that friction dominates over fluctuation as the 
source of heat. This is consistent with the fact that the 
bead is displaced from an initial stable ii -configuration to 
a final if one which is optically more confined. The small 
Qrev value quantifies the deviation from adiabaticity with 
W > A Um- This deviation stems from the mechanical de¬ 
formation of the interaction potential at both ends of the 
path which is due to an increase in the optical intensity 
as the mirror gets closer to the waist. 

In essence, our optical piston configuration enables to 
control the source of heat. This could lead to the possi¬ 
bility to reach adiabaticity with W = AUm in a simple 
way. For instance, in a pure standing wave configuration, 
both ends of the path have identical trapping stiffnesses 
leading to Q rev = 0. In this context, tailoring the optical 
landscape is particularly appealing. It leads to the pos¬ 
sibility to induce and probe all kinds of thermodynamic 
processes through the control of both heat production 
and potential energy differences. Because these quantities 
are optically determined, the level of control available is 
expected to be much smaller than k B T. 
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FIG. 3. (a) Regions of stability and bistability in the i — M pa¬ 
rameter space. Experimentally solved bistable configurations 
are displayed in open circles as i is varied. The points in the 
stable region (open squares) are extrapolated from the extreme 
bistable points according to the 20 nm piezo-actuation. The 
color map in the bistable phase codes the asymmetry of the 
bistable potential as AU/(Um(zi,£) ± Um(z 2,i)) positive in 
red, negative in blue, (b) Intensity probability densities (for 
30 s acquisition) as the mirror distance i is reduced, crossing 
over the bistability region between two stable bead dynamics. 
The six central bistable plots correspond to the six open circles 
plotted in (a). 


Appendix A: Experimental Setup 

The optical setup is sketched in FigJTJ A linearly 
polarized TEM 0 o beam from a CW diode-laser (Excel¬ 
sior Spectra-Physics, wavelength A = 785 nm, power 
45 mW) is sent into a dry objective (Nikon CFI Plan 
Fluor 60X, 0.85 NA) and focused in a water cell (deion¬ 
ized water, 80 /im thick) enclosing mono-dispersed di¬ 
electric polystyrene beads (Thermoscientific Fluoro-Max 
Red Dyed, refractive index 1.58) of radius 500 nm. The 
cell is topped by a 170 /im thick cover slip. Spherical 
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aberrations are compensated by the objective (correction 
ring set to 0.2 mm). 

The laser beam traps a single bead in the vicinity of 
a movable mirror (300 nm thick evaporated gold film on 
a glass substrate). The beam is reflected by the mirror 
and recollected by the objective. It is sent to a non¬ 
polarizing cube beamsplitter where it is equally divided. 
One arm of the intensity signal is vignetted by a pin-hole 
and recorded by a PIN photodiode (Thorlabs DetlOA). 
Amplified before numerical conversion and acquisition (NI 
PCI-6251, 16 bits resolution), this port provides intensity 
time-traces that measure the axial displacement z(t) of 
the bead inside the optical trap. The second port is sent 
to a CCD camera (Allied Guppy Pro F-031) that images 
the recollected beam spot. 

We took care to isolate optically the CW diode-laser, 
using a free-space Faraday isolator (Thorlabs IO-5-NIR- 
LP). The isolation is further improved by injecting the 
laser beam into the objective using a polarizing cube 
beamsplitter coupled to a quarter-wave plate. This pre¬ 
vents as much as possible the recollected signal to be 
send back to the injection port. This scheme implies that 
the optical landscape created between the objective lens 
and the mirror is the coherent superposition of a forward 
right-handed circularly polarized beam and a backward 
left-handed circularly polarized beam. 



FIG. 4. Schematics of the experimental setup, indicating in 
particular the secondary isolation stage at the level of the 
polarizing beam splitter (PBS) and the added quater-wave 
plate A/4. The two ports from the non-polarizing beam split¬ 
ter (NPBS) used for recording the recollected beam are also 
represented. In the chosen frame, a right-handed circularly po¬ 
larized field propagating along the z > 0 direction is described 
with cr + = (x — iy) / \/2. 


Appendix B: Solving the interaction potential 


by its Rayleigh range zr and its waist wq position fixed 
at z = 0 

E +fc (2) = Eoo T , exp (ikz - i£(z; z R )) cr+ (Bl) 
w(z;z R ) 

with 


w(z-,z R ) = w 0 y l + (B2) 

£(2; z R ) = arctan (^—^ (B3) 

and zr = ttWq/\, 712A = Ao, the optical wavelength in 
water (refractive index 77,2) and ttwq = A/NA. 

The reflected field, counter-propagating k —k with 
respect to the incident beam, is multiplied by a reflectivity 
coefficient r(A). Its waist position is the mirror-image of 
the incident waist position. But before being reflected, the 
incident beam is intercepted by the bead because the bead 
diameter is larger than the incident waist wq. Considering 
that the refractive index of the bead (polystyrene) is 
different from that of the fluid (water), the incident beam 
transmitted through the bead is magnified, the bead 
considered to act as a lens-doublet. We account of this 
effect by introducing an effective magnification parameter 
M on the reflected beam itself. This leads to changing 
the Rayleigh range zr — >> M 2 zr and waist wq —>> Mwq 
of the reflected beam with respect to the incident beam. 
The reflected beam, left-handed cr_ circularly polarized, 
is then expressed as: 


E 


Mw 0 


M ( z ^£) — ^00 / , ,2 \ 

w(Z]M z Zr) 


exp [—ik(z — 2£) + i£(z — 2£; M 2 zr )) <r_. (B4) 


The coherent superposition of the incident and the 
magnified reflected Gaussian beams determines the optical 
landscape of the problem. The corresponding optical 
intensity writes as 


= \E+ k (z) + r[X]-E- M k (z-2£)\ 


— F 2 ^0 

00 w 2 (z;z r ) 


+ p 2 e ( 2 


00 


''0 


+ 


2pEl 0 Mwl 


w 2 {z-2i-M 2 z R ) 

X 


(B5) 


w(z; z R )w(z - 2 £; M 2 z R ) 

cos (2k(z - £) - £(z; z R )~ 

£(z ~ 2^; M 2 z r ) + V>) • 


1. Optical Fields 

The optical fields E to t created inside the optical piston 
(i.e. between the objective lens and the mirror) is given by 
the coherent superposition of an incident and a reflected 
Gaussian beams. The incident beam is, as explained 
above, right-handed cr + circularly polarized and described 


with an interference term between the two beams. In the 
vicinity of the waist, it shows that the modulations of 
the optical landscape (that will eventually correspond to 
the local potential barriers, as discussed below) are deter¬ 
mined from an harmonic term. This immediately stresses 
that a standard 4 th order polynomial description of the 
barrier -as done when resorting to a typical Duffing model- 
is not appropriate for our optical piston configuration. 
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2. Conservative optical force 

In the dipolar regime, the bead is characterized by 
an electric polarizability a. Neglecting any source of 
dissipation within the bead, i.e. assuming that lm[a] ^ 0, 
the gradient force is the only force exerted on the bead 

F op t = -d z U M {z,l). (B6) 

It directly derives from the interaction potential energy 
determined from the optical landscape intensity as 

Um(zJ) = -A • — I op t{d,z) (B7) 

2 e 0 n 2 c 

with a = 47reori2a 3 [( 711 / 77 , 2 ) 2 — l]/[(ni/ri 2 ) 2 + 2], n\ the 
refractive index of the bead, 712 the refractive index of 
the fluid, a the radius of the bead and eo the vacuum 
permittivity. 

Despite the fact that is is related to the bead geometry, 
the M parameter value for the magnification by the bead 
is not necessarily a constant of the model as it depends 
on the width of the beam intercepted by the bead. The A 
parameter quantitatively corrects the value of the poten¬ 
tial calculated in our Rayleigh-based model due the finite 
size of the bead that should be accounted for in a more 
realistic description of the optical interaction. In fact, as 
given in the main text, the interaction potential turns 
out to be smaller by approximately 3 orders of magnitude 
for a 1 /im bead. This value is in good agreement with 
existing evaluations for the size effect |13j . We stress that 
because A is a parameter of the model that characterizes 
the intensity of the coupling of the bead with the optical 
intensity, it is kept fixed once determined for a given bead. 


3. System Resolution 


Monitoring the instantaneous motion I(z(t)) of the 
bead in the bistable phase, mean residency times and 
stiffnesses are extracted from experimental intensity 
time traces for each well. Kramers rate equations provide 
the following system fed with the extracted values: 


d 2 z U M (z,£)\ z=Zl 

d 2 z U M (z,£)\ z=Z2 

< \A^l\/^b ( UM(zb,£) ~ Um(zi,£) \ 

2ttt] GXP \ k B T ) 

\/^2\/^b ( Um[z^) — Um[z2J) \ 

„ 27T7? exp \ k B T ) 



(B8) 

K2 

(B9) 

1 

n 

(BIO) 

l 

T2 

(Bll) 


where Eqs. (BIO) and (Bll) are escape equations given 
by Kramers theory in the over-damped regime for each 
of the metastable states of the bistable phase. Resolution 
of this system gives access to the remaining unknown 
quantities of our model: the mirror waist distance £, 
the bead magnification effect M and the size correction 


parameter A related to the coupling between the light 
field and the finite size bead. 

The energy difference between the two metastable wells 
can be derived from the extracted residency times and 


stiffnesses as A U = &bT In • Taking the ratio of 

equations |B10| and |B11| then removes the dependency of 
the system on the properties of the barrier but requires 
the knowledge of A. Fixing A and using the simplex 
algorithm, the external variable i (piston length) and 
the magnification parameter M are determined. This 
therefore gives the entire potential Um{z,£) for a given A. 
Iterating this resolution over A until the rate equations 
are verified provides the triplet (£, M, A) that best solves 
the whole system. Once determined in one ^-configuration 
(say the symmetric bistable configuration), the parame¬ 
ter A is then kept constant when solving other bistable 
configurations (for instance varying £). 


Appendix C: Axial Displacement 

The bistable behavior of the bead can be monitored on 
the CCD camera because the motion of the bead along 
the optical axis and across the bistability barrier changes 
the Gaussian envelope of the reflected beam and thus the 
diffraction pattern of the recollected beam imaged on the 
camera, as shown in Fig Thus, by the sole measure¬ 
ment of the recollected intensity, one can access part of 
the bead dynamics. Nevertheless, the low acquisition rate 
of the CCD camera is a strong limitation for analyzing 
precisely the stochastic motion of the bead. 


(a) 

(b) 

<• 

» 


FIG. 5. Recollected spots imaged on the CCD camera in the 
case of a bead trapped in a bistable configuration, with same 
exposition time. Panel (a) shows the bead in its most distance 
position from the mirror, i.e. the position behind the waist 
(z < 0). Panel (b) shows the bead in front of the waist, thus 
closer to the mirror. Both images have been recorded with 
the same exposure time. The recollected intensity is higher 
near the mirror (image (b)) than away from it. The central 
area of the recollected spot is indicated by the superimposed 
circles, with a 44 pixel-diameter for (a) and 50 for (b). 

To do so, we rather measure the recollected intensity 
using a PIN photodiode that grants a better sensitivity 
and a high acquisition rate. We set this rate for our 
experiments at 2 18 = 262144 Hz using low noise preampli- 
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fiers (SR560). The PIN signal for intensity measurements 
was recorded in AC mode thus filtered through a 0.3 Hz 
high-pass filter at 6 dB/oct to remove the continuous com¬ 
ponent of the signal. High-pass filtering poses no issue 
since we focus on signal fluctuations while it allows, after 
amplification, to span the signal of interest over the whole 
acquisition card input range. For the slow varying signals, 
the exponential decrease coming from the high-pass filter 
is later compensated numerically to within a constant. A 
low-pass filter at 100 kHz at 6 dB/oct was also used to 
avoid aliasing. 

Time traces measured by the PIN photodiode clearly 
reveal the bistable dynamics of the bead as shown in Fig [6] 
The intensity time trace is distributed around two mean 
values from which we extract and concatenate the dynam¬ 
ics associated with each well of the bistable potential, as 
shown on Fig [ 7 ] The full signal dynamics clearly appears 
stationary: the signal stays centered around a constant 
value (near 0) and the ratios between the mean residency 
times for each mean value is time-independent. Simi¬ 
larly, the two concatenated time traces reveal stationarity 
properties. 



0 10 20 30 40 50 60 


Time (seconds) 

FIG. 6. Time trace associated with a bistable motion of the 
bead as recorded by the PIN photodiode on the reflected beam 
intensity. The acquisition time is set to 60 s, at a rate of 262 
kHz. 

In the experiment, we spatially filtered the recollected 
beam before the PIN photodiode through a pinhole (typi¬ 
cal aperture of about 1 mm 2 ). This filtering was useful 
since it actually enhanced the separation between the two 
average intensities in the bistable configurations. 


Appendix D: Power Spectral Density 

Assuming that the recorded stationary intensity time 
traces are ergodic, meaning that a time trace is indepen¬ 
dent of the initial position of the bead and leads to the 
same distribution for different realizations with identical 
parameters, one can study the bead dynamics through 
its power spectral density (PSD) [21 j. This standard 


0.3 



0 5 10 15 20 25 30 35 

Time (seconds) 

FIG. 7. Concatenated time traces associated with each of the 
mean values measured in the time trace of Fig. [6] 




FIG. 8. PSD taken on each of the concatenated time traces 
given in Fig. [ 7 ] Upper curve is the PSD for the I(t) > 0 
concatenated time trace and the bottom one the PSD of the 
concatenated time trace when the bead is in the well further 
away from the mirror (corresponding to I(t) < 0). The black 
solid line is the Lorentzian fit of the data and the vertical line 
gives the roll-off frequency of each of the quasi-harmonic trap 
associated with each locally stable positions. 


approach has the advantage of being straightforward to 
apply on rapidly fluctuating data. 

It clearly appears that the PSDs associated with each of 
the concatenated time traces of Fig. [ 7 ] follow the typical 
Lorentzian shape of a Brownian motion performed in 






















































FIG. 9. PSD taken directly from the full time trace shown in 
Fig. [6] for an acquisition time of a minute. The bottom PSD 
(dark) corresponds to the experimental laser noise. 


an harmonic trap in the over-damped regime, as seen in 
Fig. H This implies that the two wells of the bistable 
potential, separated by the activation barrier, are quasi¬ 
harmonic, with stiffnesses ftp 2 = 27r 77 /p 2 that can be 
determined directly from the fluid drag r] and the so- 
called roll-off frequency /p 2 of the trap measured on each 
PSD [14]. In the full signal PSD shown in Fig. [9j these 
quasi-harmonic traps are seen through the Lorentzian 
fit at high-frequencies. But the spectrum within these 
local wells does not exhaust the bistable dynamics. Low 
frequencies indeed reveal a strong increase in the power 
spectrum which is due to the activation process of the 
bead over the bistable barrier, occurring on a typical ~ 1 
Hz regime. In other words, the crucial separation of time 
scales discussed in the main text is readily observed on 
the PSD associated with the bistable motion of the bead. 


Appendix E: Experimental Uncertainties 

We assume that the distribution of the measured res¬ 
idency times are Poissonian, implying that their mean 
values rp 2 equal their variances oy . Measuring N = 25 
number of back and forth activations of the bead through 
the barrier (over an acquisition time of T = 30 s) for the 
bistable configurations described in the main text, leads 
to an experimental uncertainty in the Tp 2 determination 
of St p 2 = g Ti 2 / y/N = 0.2 0 > 12 . Because the signal is 
stationary, t\ + T 2 = ^ and therefore 72 = ^ — t\ . 

In addition, by taking the expression of the perpendic¬ 
ular viscosity [22] and its derivative, and neglecting the 
systematic error in the region where the bead evolves, the 
uncertainty in the change of viscosity along the displace¬ 


ment of the bead can be estimated at a 3% level over 300 
nm bed displacement. To this uncertainty, 1 % is added 
from the extraction of the roll-off frequency at the level 
of the PSD (fit uncertainty). This thus leads to a global 
stiffness uncertainty of about Sk p 2 = 0.04 ftp 2 . 

The 6 resolved M values have mean and sample de¬ 
viation respectively of M = 1.508 and gm — 0.071. 
The uncertainty therefore is SM = gm/V$ = 0.029 and 
M = M dz SM = 1.508 zb 0.029 which is an uncertainty of 
4% of the value. 

These uncertainties are propagated [23J to determine 
the uncertainty on A U used as an input (through fti ? 2 
and T 12 ) in our system solver. Propagations of 5t\ and 
Sk i ? 2 in A U in bistable configurations where t\ ~ r 2 lead 
to 



dAU 

dr\ 


(M) 2 + 


(dA U\ 

\~d^J 


= <°- 2 > 2 ( 1 + f +2 b + 
~ 0.16. 


2 (fc ' ) 2 + (^r ) !<&2)2 

(°?P(T ) 2 


Because the propagation of uncertainties is logarithmic 
for SAU, resolution of energy differences between the two 
well is lower than half a k^T despite the uncertainty on 
average lifetimes. 

Taking M = M (given the small SM value) allows us 
computing the sensitivity of AU as a function of t. The 
computation yields dAU/dt = 0.07 k^T • nm -1 . Combin¬ 
ing this sensitivity with the uncertainty SAU of 0.4 k^T 
gives an uncertainty on the waist-mirror position St of 6 
nm only. This rather high spatial resolution is an inter¬ 
esting by-product of our approach. 

The reversible heat measured with our method on our 
experimental configuration through the cross-over path 
(see main text) is computed from trap stiffnesses which 
depend on the waist-mirror distance t. The measured 
heat uncertainty SQ rev produced along the path is thus 
estimated from the determination of Si. The trap stiff¬ 
ness of the stable positions (around 2 = 0) is computed 
for an incremental displacement of i of =b St. A worse- 
case scenario is then followed, taking the highest dif¬ 
ferences in trap stiffnesses between k(£) and + Si). 
The heat uncertainty can then be computed as SQ rev = 


k^T In y/K,{jt)/K(t + St) ~ 2.2 x 10 23 J with ft(£) = 

6.47 x 10“ 6 N • m- 1 and + St) = 6.40 x 10“ 6 N • mW 
Similarly, the worse-case uncertainty for the potential 
energy SAUm = Um (+ St) — Um( 0,^) is determined 
and is about 6 x 10 -23 J. 
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